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Abstract 

Chiral-odd structure function hi is expected to be measured in polarized 
Drell-Yan process. We calculate two-loop anomalous dimensions for hi in the 
minimal subtraction scheme. Dimensional regularization and Feynman gauge 
are used for calculating the two-loop contributions. Our results are important 

in studying dependence of hi. 
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Internal spin structure of the proton has been a popular topic since the discovery of a 
serious proton-spin problem at high energies. In spite of a naive quark-model prediction, it 
was suggested that the major part of proton spin is not carried by quarks by measuring the 
structure function gi. The discovery inspired many theoretical and experimental studies on 
the spin structure. In order to understand how the proton spin consists of quark and gluon 
spins, it is important to learn about sea-quark and gluon polarizations. As an experimental 
project to investigate these issues, the RHIC (Relativistic Heavy Ion Collider) Spin project 
was proposed. 

In this project, the internal spin structure will be investigated by polarized proton-proton 
reactions. Another virtue of this project is that a chiral-odd structure function hi could be 
measured in the Drell-Yan process with transversely polarized protons M. It could be 
also measured in semi-inclusive electron scattering . It provides us another important test 
of proton spin structure. In particular, the structure function hi probes different aspects 
of spin. Furthermore, a nice point is that it has a clear parton-model interpretation. It 
measures the transversity distribution, namely the probability to hnd a quark with spin 
polarized along the transverse spin of a polarized proton minus the probability to hnd it 
polarized oppositely. Its gross properties have been studied recently in particular by Jaffe 
and Ji |0|. Here, we would like to focus on dependence. Leading-order evolution was 
already studied p|. However, it is our standard to use next-to-leading-order (NLO) results 
in obtaining unpolarized and polarized parton distributions. Recently, the NLO evolution 
equations are established for the structure function gi |0. Therefore, it is important to 
investigate the NLO evolution of hi for future analyses of polarized Drell-Yan data in 
connection with the hi structure function. Our research purpose is to calculate two-loop 
anomalous dimensions, which determine the NLO evolution of hi. 

In order to study hi in perturbative Quantum Chromodynamics (QCD), we need to 
introduce a set of local operators 0 

i 75 — trace terms , n = l,2,... . (1) 
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Here, Sn symmetrizes the Lorentz indices /ii, /i^, and iD^ = id^ + gt°‘A^ is the covariant 
derivative. We calculate two-loop anomalous dimensions of these twist-two operators. Be¬ 
cause gluon held cannot contribute to the anomalous dimensions of the chiral-odd operators 
even in the NLO, our calculation procedure becomes a rather simple one which is similar to 
the nonsinglet formalism for example by Floratos, Ross, and Sachrajda |Q. Renormaliza¬ 
tion of the '^7^1 (1 ± 75)zZi)^2 ... ^^ypg operators was studied in Ref. |P], but the same 

method can be applied to the operators in Eq. i). First, Feynman rules at the operator 
vertices should be provided. In order to satisfy the symmetrization condition and to remove 
the trace terms, the tensor with the constraint A^ = 0 is usually multiplied. 

However, the operators in Eq. ([^) are associated with one-more Lorentz index z/, so that it 
is convenient to introduce another vector with the constraint ■ A=0. Then, Feynman 
rules with zero, one, and two gluon vertices become those in Fig. |I|. 

Next, we give a brief outline how anomalous dimensions are calculated |^. The bare 
operator is dehned by = V'sQbV'b = with the renormalized one Then, the 

renormalization constant Zq is given by Zq = Zp/Zg, where Zp and Zg are wave-function 
and operator renormalization constants. Anomalous dimensions for O"', Q"', and 'ijj are 
given by these renormalization constants as = gid{\nZo^)/dgi, = g.d(\nZgn)/dfi, 
and 7 ^ = {fi/2) d{lnZp)/dfi. Using the relation among the renormalization constants, we 
express the anomalous dimension 7 ^^ in terms of and 7 ^ as 7 ^^ = 27 ^ — 7 ^^. In 
dimensional regularization with the dimension d = 4 — e, Zgn can be obtained in the form 
= 1 + S^i As it is shown in Ref. 0, the anomalous dimension is related to 

the hrst coefficient Z^^^ by 7 ^^ = —g^oZ^ /dg^. In the two-loop calculation, it becomes 
7 q „ = — 2 Zi . The anomalous dimension 7^„ is simply given by minus twice the 1/e 
coefficient in The same equation can be also applied to the anomalous dimension 7 ^. 
In this way, we try to hnd 1/e singularities in the renormalization constants for calculating 
the anomalous dimensions. 

Two-loop contributions are calculated by using the dimensional regularization with the 
Feynman gauge. Here, minimal subtraction (MS) scheme is used. Possible Feynman di- 
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agrams contributing to the one-particle irreducible matrix elements are shown in Fig. |^. 
Furthermore, there are two-loop diagrams contributing to the quark self-energy as shown in 
Fig. 1^. However, the quark-held renormalization is already discussed in Ref. P, so that we 
do not have to repeat the same calculations. In order to explain our calculation procedure 
for the operator part, we pick a simple case. In Fig. Lorentz and color indices are shown 
with momenta for the diagram (d) in Fig. |^. The diagram is calculated with a Feynman 
rule in Fig. |l| as 


Ih = 


d^ki d^k2 . ^ iJ^2 -iguphc ■ ,c p 

ptp(2tF^ ^ IF (FTVfeF^ F 

n—1 

X ^ 75^^ (A ■ py-^ gt^/V (A ■ . 

i=i 


( 2 ) 


Calculating the above integral, we obtain 




1 r(e/2)F(6) 

( 47 r )4 F(1 + e/2) 


R(2-e/2,1-e/2) 


n—1 


X H B{j + l-e,l-e/2) 
i=i 


—p 

dvr 


(3) 


where F(x) and B{x,y) are gamma and beta functions, and Cp is given by Cp = {N^ — 
l)/{2Nc) with the number of color N^. Subtracting a one-loop counter term from the above 
equation and evaluating singular terms, we have 




n—1 


(dvr) 


clU- 


J + 1 L e 


- {1 - Si{j + 1 )} 


(4) 


where S'i(n) = 1/^- From the 1/e singularity in the above equation, we obtain the dia¬ 

gram (d) contribution to the anomalous dimension as 7 ^^'^^ = ^9^ Cp [ Gi (n) -S'i(n)]/(d 7 r)^ 
with Gi{n) = I]j=i(l/j) Si=i(!/*)• The factor of two is included by considering a similar 
diagram with gluons attached to the initial quark line. We hnd that the obtained anomalous 
dimension is exactly the same with the one in Ref. . The reason is the following. Because 
the operator vertex part 75 ^^/! can be separated from the ki and k 2 integrals in Eq. (^, the 
integrals are independent of the operator form. The vertex in Ref. is simply replaced 
by the present one 75 ^^^. Therefore, the renormalization calculations are the same in both 
cases, so that the obtained anomalous dimensions are exactly the same. 


d 












We notice that several diagrams in Fig. ^ have the same property in the sense that the 
operator part can be separated from the integrals. These are cases with gluon propagators 
attached only to the hnal-quark or initial-quark line. In addition to the diagram (d), these 
are diagrams (g), (h), (1), (m), (q), and (r). We do not have to repeat the same calcula¬ 
tions. Furthermore, there are certain diagrams which do not contribute to the anomalous 
dimension. Calculating diagrams (a), (b), (c), and (p), we hnd that there exist no 1/e 
type singularity. Therefore, there is no contribution to the anomalous dimension from these 
diagrams. From the above discussions, the problem reduces to calculations of remaining 
diagrams (e), (f), (i), (j), (k), (n), and (o). 

We need to calculate integrals which appear in evaluating these diagrams. Fortunately, 
many of them are already listed in Appendix of Ref. [§]. All divergent subintegrals are 
renormalized by subtracting the pole part of the subintegral. However, special care has to 
be taken in some diagrams, for example the one in Fig. ^(i). Contracting 7 matrices in 
the numerator, we hnd a term proportional to e. In this case, the 1/e contributions come 
from the double-pole terms, in which it is known physically that there is no counterterm 
contribution [Q. Therefore, when we write the double integral results, it is more convenient 
to express the 1 /e^ terms without subtracting counterterms even though the 1 /e terms 
are calculated with the counterterms. This is actually done in Appendix of Ref. [§ for 
writing the integral results. We also express extra integrals in this way in Appendix B. 
Once we have all the integral information, it is straightforward calculation for obtaining the 
anomalous dimensions. Because calculation procedure is essentially the same with the one 
for the diagram (d) and the details are already described in Ref. 0, each evaluation is not 
explained. We simply list our results in Appendix A. 

We summarize our study. Two-loop anomalous dimensions for the chiral-odd structure 
function hi are calculated in the MS scheme. Dimensional regularization and Feynman 
gauge are used for calculating the possible diagrams. Our results are important for studying 
the details of dependence in hi, which is expected to be measured at RHIC. 

Note added: After we submitted this paper (|hep-ph/9706420|) , two preprints appeared 
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in the Los Alamos archive on the same topic. Our anomalous dimensions agree not only 
with the Vogelsang’s result (|hep-ph/9706511| ) but also with the Hayashigaki, Kanazawa, and 
Koike’s (|hep-ph/9707208D although the integral expressions are different in some diagrams. 


Appendix A: List of anomalous dimensions 

Contributions from the diagrams in Figs. @ and ^to the anomalous dimension are listed. 
In the following equations, Hk is dehned by Hk = b = 0)^—Hk{n—l, b = 1)4^], 

and DP (SP) denotes double-pole (simple-pole) part of the integral. The integrals Hk are 
listed in Appendix B, and A are the integrals given in Appendix of Ref. |^. The function 
F is defined by F{a, c) = c! F{a + J,c— J)/[J\ (c — J)!]. Some results are 

exactly the same with those in Ref. |§] as explained in the text: 

^( 2 rf) ^( 29 ) i2h) (21) (2m) {2q) (2r) 

7i“’, 7^', 71“’, 7^’ given ill Ref. [8] . (5) 

Other anomalous dimensions are calculated as 


H2a) ^ (2b) ^ (2c) ^ (2p) ^ g 

in In In In ^ 


( 6 ) 


I n 


9 ^ r T - I r r 


(7) 


^ = + 9 " - CfCg/2) 4 (A ■ 


4 - 6 H^^{n) + 4 H^^{n) - 4 H^^{n) 


2 If ^(0, n - 1) - (0, n - 1) + 21,^(0, n - 1) - 1,^(0, n - 1) + 3 1^/(0, n - 1) 


tDPi 


tDP 


DP/ 


tDP/ 


+ (A-p) Jf^(0,n) -21^^{l,n-l) - n) + n - 1) +n) 

n - 1) - /f/(0, n) + 2 n - 1) I + 4 Hl^{n) - 21^^{n - 1, 0) 

+ (A-p)-i| -2Jf^(0,n) + 2Jf^(0,n) -4/f^(l,n- 1) + 2/3^^(n,0) 

+2 - 1,1) + 4 I^in) - 2 /f/(0, n) + 2 n - 1) + 2 4^(0, n) 

-2 Jf3^(l, n - 1) - 2 /f/(0, n) + 2 /f/(l, n - 1) | 

+ (A ■ p)-2| 4 + 1, 0) - 8 1) + 2 /f^(0, n + 1) - 2 J|^(0, n + 1) 
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+4n) - 2 I^^{n + 1, 0) - 4 jf^(n + 1) - 4 I^^{n + 1) + (-l)M 7f^(n - 1, 2) 

+2 /f/(0, n + 1) + 2 /f/(2, n - 1) - 4 n) - 2 4^(0, n + 1) 

+4 n) - 2 I^f{2, n - 1) + 2 n + 1) - 4 n) } 

^ “ C'fC'g/2) 8 , (8) 

^(2i) = + / C'^C'^ (A ■ 8 -^(n) - 4# 2 ^^(n) + 8 ^ 3 ^^(n) + Jf ^(0, n - 1) 

+/f ^(n - 1, 0) - 4n - 1) - 5 - 1, 0) + (A ■ p)-^| - 3 Jf-^(0, n) 

+4 1^(1, n - 1) - Jf ^(n - 1,1) + 3 1^(0, n) - 4 Jf ^(1, n - 1) + -1,1)} 

-8 i72^^(n) + 2 /f-^(0, n - 1) + 2 /f^(n - 1, 0) + (A ■ p)-^| 6 /f^(0, n) 

-6 /f^(n - 1,1) + 4 /f^(l, n - 1) + 2 n) + 6 I^^{n - 1,1) 

-8 /2''"’(1, n - 1) - 8 0) + 4 li^{n - 1,1) - 8 I^^{n) } 

+ (A ■ p)-2| (-1)”+46 1) + 4 /f-^(n - 1, 2) - 8 /f^(l, n) - 4 I^^{n - 1, 2) 

+8 / 2 ^^( 1 , n) - 8 /f^(n + 1) + 8 + 1) - 8 I^^{n -1,2) + 16 1) | 

= 0 , (9) 

^(2fc) = +g^ (6-2 - CfCg/2) 2 (A ■ p)-(-i) [ (-1)" 8 ^(n) - - 1, 0) + 2 I^/{n - 1, 0) 

- 1, 0) + 3 - 1, 0) + (A ■ p)"^| I^^{n - 1,1) - /iT(n, 0) 

+Iu^{n -1,1) + /iT(n, 0) - - 1,1) - I^/{n - 1,1) | 

+8 {1 - (-l)'^} {Hrin) - + (A ■ p)-^I^AH^^{n - 1,1) -AHf^{n, 0) 

-4 H^^in - 1,1) + 2 Ii^{n - 1,1) - 2 0) - 2 7f^(n, 0) + 2 /f/(n - 1,1) 

-2 /f/(n, 0) - 2 I^fin - 1,1) + 2 ^^(n, 0) - 2 I^f{n - 1,1) | 

+ (A ■ p)-2| (-1)^ 4 H^^{n + 1, 0) + 8 1) - 4 H^^{n + 1, 0) 

-4 iy7^^(n - 1, 2) + 4 i77^^(n - 1, 2) + 2 I^^{n + 1, 0) - 4 I^^{n - 1, 2) - 8 7|^(n, 1) 
+47|^(n-l,2)| 

= - tTs (+-8 {1 - (-1)"} ^Trn • <1°) 

(47r)^ n(n + l) 
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n —1 p P 

7 '“’ = + 4g" CfCg (A . p)-<“-» - 1 - g) + (A . p)-> - (n - 1 - j, j + 1 ) 

i=i ^ 


+A U n- i) + l 2 {n-l- j,j + 1) - /a (j, n-j 

,4 r X 




( 47 r) 


C'fC'g4 


{4^i(n) - _ ^ 2 (n)} + S,in) - 2 S{n) 


( 11 ) 


n—1 r 


7 ‘"°> = + 9 * (C| - Cp(7g/2) 8 (A ■ p)-<"-i) 


i=i 


AT (^ -1 - j, j) + AT i -1) 




+/f/’(n - 1 - j, j) + /f/(n - jj - 1 ) + (A ■ p)-^| + /T(p -l-jj + 1) 
-/T(p - j, j) - AT -1 - J, J + 1 ) + 2 AT- A j) - AT- J + 1 , J - 1 ) 
+AT(^ -1 - j, j + 1) - 2 AT j) + AT- j + 1, j -1) 


- A 4 (p -1 - J, J +1) + A 4 - j) i) 
-4 ^ ^ 1 




(Att) 


{Cl - CfCg/2) 8 2 ^i(n) I - - 2 S'^{n/2) + 2 ^ 2 (p) } 

-2S^{n) - S'^{n/2) + 4^(n) + SS{n) 


( 12 ) 


In these equations, the color constants are Cg = Nc^ Cp = {N^: — l)/(2Ac), and Tr = Nf/2, 
where A,, is the number of color and Aj is the number of flavor. The functions Si, S 2 , and 
others are dehned by S'^p) = YJk=i 1/A, S''(n/ 2 ) = [{!+(—1)"'}S'j(p/2) + {1 — (—1)”}S'A(p— 
l)/2)]/2, S{n) = Z]fc=i(—l)^*S'i(fc)/fc^, and S{n) = YCk=i Si{k)/k‘^. We have to be careful in 
handling the 13 ( 0 , 0 ), / 8 ( 0 ,c), and ^ 7 ( 0 , 0 ) terms; however, the above results are valid also 
for n = 1. From the above expressions together with other diagram results including the 
quark-held renormalization part in Refs. Q] and [Q, we obtain the anomalous dimension for 
the structure function hi. 


^ total 
In 


= Cl 


32^i(n){^2(n) - T(T2) } + 24R2(T - 8S{i{n/2) + 64F(n) 
ol-(-l)’ 


- 3 


+ CpTp 


160 32 4 

-^>5i(n) + —S2{n) + - 


+ CfCg 


n{n -|- 1 ) 

p:o^ 00 

-^Si{n) - -^S 2 {n) + AS',{n/2) - 32 TT 
y o 


lQSi{n){2S2{n) - ^AT2) } + 4 


l-(-l)^ 

n{n + 1 ) 


17 

y 


(13) 


where the factor 5 f^/( 47 r)^ is taken out. 
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Appendix B: Integrals 

We list integrals H^. The simple-pole terms are written with connter terms and the 
donble-pole terms are withont them. Finite contribntions are neglected. Becanse the H 
fnnctions are nsed for Ffi-s in Appendix A, we express their resnlts in the H form: 

H ( h) = f f _ (A ■ kiY {Q ■ _ 

^ ’ J { 2 TrY J { 2 TrY {kiY {ki — pY {k 2 — kiY {k 2 — ki + pY ’ 

^i(«) => (14) 

where 5i(e, a) = —4/e^ — (2/e)/(a -|- 1). The following integrals are given by the above 
expression with some replacements: 

H 2 '. {k 2 — kiY —^ ^2 o) = —2/e^ -|- (l/e)[—l/(a -l- 1) -l- 1/a — 2], 

H^\ {k 2 - ki +pY kl and ^i(e,a) = -2/e^ (l/e)[S'i(a) - l/(a 1) - 2], 

H^\ (fci - pY{k2 -ki+ pY {k2 - pY{p + ki- k2Y 

and 5i(e, a) = -{2/eY/a + (l/e)[l/(a + 1) -l/a- l/a?], 

H^\ ^i/[(fci - pYih - ki +pY] h/[ik 2 - p)‘^{p + ki - k 2 Y]- is expressed by 

Eq. (U) with ^i(e, a) = -|-(l/e)/a. 

Different integral types are 


HYa,c) = 


dk^ki f d‘^k2 


(A ■ kiY (A ■ k2Yk2 ■ p 


{2'kY j (27r)‘^ (ykiY {ki — ^2 
1 (-l)"a!c! 


{At^Y (a -|- c -|- 1)! (a -|- c) 


2 kl (^2 - ki +pY 

(A.pr+'[ -i+l{-S,(a + c) + Si(a) 

+Si(c) -i—--1 + 1 

d H- c 1 Qj 


(15) 


r (Yki r (Yk 2 ■ kiY ■ k 2 Y ki ■ p 

^ J {2'kY J {2kY {kiY {k2 - pY (kl - k2Y {p + kl - k2Y 

1 /A Na+c^ (Q + c-^-l)!(^+l)!a!c! r 1 
(4^)4l -P) J)!(a + c)!(a +J+1)! [ 

— 25 * 1(0 -|-c-|-l)-|-S'i(a-|-c — J — l)-|-S'i(a-|-t7 -|-1)|’ 


(16) 
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FIGURES 



Yj jX (A-k)" ■ ‘ X'j %t!L0. (A-k,) J‘ g f A^‘i (A-k2)"■ ‘ ■ -i 



g2Y5IX^2^A^‘A'' £ £ (Aki)‘-‘ [ t'>t“{A(k2+q2)>'-‘ 

+ t‘‘t'’{A(kj-q2))J-‘] (Ak2)"'^-J 

FIG. 1 . Feynman rules for zero, one, and two gluon vertices. 
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FIG. 2. Contributions to the two-loop anomalous dimension. 
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(a) (b) (c) 


vL/V 

(d) (e) 


FIG. 3. Two-loop contributions to the quark self-energy. 


o,d 



FIG. 4. Lorentz and color indices are shown with momenta for the diagram in 


Fig. 1(d). 
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